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TRANSONIC FLOW PAST AN AIRFOIL 
WITH CONDENSATION 


B. Schmidt 


SUMMARY 


/ 194 * 


In connection with investigations conducted to determine the influence of water 
vapor on experiments in wind tunnels, the question arose as to what changes due to 
vapor condensation might be expected in airfoil measurements. 

Density measurements on circular-arc airfoils aided by an interferometer in 
choked tunnels with parallel walls show that increasing humidity produces increas- 
ing changes in the flow field. The flow becomes nonstationary at high humidity. 

However, the influence of the condensation is only felt at the airfoil, inasmuch 
as the shock bounding the local supersonic region moves upstream with increasing 
humidity while its intensity decreases. The density distribution upstream of the 
shock remains unchanged. Even if the flow becomes nonstationary in the vicinity of 
the airfoil, no changes occur at the airfoil - the base of the shock included. This 
is a somewhat surprising and unexpected result. 


INDEX OF NOTATIONS USED 


P — p<n n 

e„ = -y — « - 2 (p, 

Qx> ti'o» 

C P — T V. * 

d 

H 


- speed of sound 

- pressure coefficient 

- standardized pressure coefficient 

- maximum airfoii thickness 
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- tunnel height (from airfoil centerline to tunncJ. 
wall) 


^Numbers in margin indicate pagination in foreign text. 


L 

M 

P 

u 

X 

x [sic] 


4 > 

4> x - 3<j>/3x - u/c* - 1 £ M* - 1 

X 

P 

T = d/L 
£ ■ x/L 
Indices 


- length of airfoil 

- Mach number 

- pressure 

- velocity in the x-direction 

- airfoil abscissa (station) 

- absolute humidity of the air (in kg of water per 
kg of dry air) 

- perturbation potential 

- perturbation velocity in the x-directiou 

- ratio of specific heats 

- density 

- normalized airfoil thickness 
normalized airfoil coordinate (station) 


0 - rest conditions 

<» - free-stream (unperturbed) conditions 
* - referred to conditions at M = 1 




I. INTRODUCTION 


Based on the research done at this institute on the influence of air humidity 
on the flow field in a Laval nozzle [Refs. 1, 2], the question arose whether simi- 
lar Influences of humidity occur in the flow about wings. Observations of con- 
densation phenomena for wings have been made [Ref. 3] but do not shed any light on 
the questions raised here. 

Simple preliminary investigations using Schlieren photography clearly show /196 

(Figures 1, 2) that, for a flow about a circular-arc airfoil represented by a half- 
model, the effect of air humidity is to change the shock configuration at the wing. 
Subsequent studies yielded quantitative results on the details and are described 
below. In particular, high air humidity can lead to nonstationary flow phenomena 
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in the vicinity of the wing. This, in turn, might lead to nonstationary forces act- 
ing on the wing, which is important in practical applications. 

Construction of the Test Facility 

The measurements and photographs were made in a wind tunnel at the Institute for 
Flow Studies and Flow Mechanisms of the Technical University at Karlsruhe, designed 
for intermittent operation (diagrams in Figures 3, 4).* The pressure drop required 
for the operating of the tunnel was obtained by pumping down a bank of tanks with a 
total volume of 34 m^. Air of a given humidity was stored in polyethylene balloons 
with a volume of about 18 m^ located upstream from the intake tube. The test sec- 
tion proper had parallel walls and measured between 35 x 100 mm and 35 x 150 mm and 
was within the field of view of the optics. The circular-arc airfoil in the test 
section was either a half-model attached to the upper tunnel wall or a symmetrical 
wing at zero angle of attack located in the center and glued to the observation 
ports. Upstream from the airfoil, the boundary layer could be sucked off through 

/ 19 7 

seven slits 1 x 40 mm each, which terminated 10 mm in front of the wing. The amount 
of suction could be held constant during the observation period and was chosen so 
that all but about 0.25 mm of the boundary layer was removed. A variable nozzle 
installed about 0.7 meter downstream from the test section was used to control the 
Mach number of the incident flow. The Mach number for the wing was varied from sub- 
critical values up to the choking of the tunnel. For t = 0.1, H/L ■= 2, and 
H = 150 mm, the corresponding range was from Mcccj-iticai £0.78 and M oochoke flow £ 

0,'83. 

A manually operated ball valve with a nominal size of 200 mm served as the shut- 
off device for the tanks and could be opened or closed in 0.2 second. The follow- 
ing parameters were measured; Static pressure using a pitot tube and static tube; 

*1 am grateful to Mr. H. Lehmann for his help in conducting the research and evalu- 
ating the photographs. 
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static pressure at the side walls along the axis of the tunnel and at 15 to 20 mm 
upstream from the edge of the field of view; static pressure at the half-model 
attached to the wall for values of £ * 0.25, 0.5, and 0.75; density distribution in 
the field of view using a Mach-Zehnder interferometer; and air humidity of the sup- 
ply air in the balloons by using an aspiration psychrometer . 

All pressures were measured by mercury U-tube manometers, whose connecting tubes 
to the test section were opened for only about 10 seconds during the stationary flow 
phase. The Schlieren and interferometer photographs were all made using spark pho- 
tography with an exposure time of about 10“^ seconds. At the times when nonstation- 
ary processes were expected, slow-motion films with a frame frequency of 10^ frames 
per second were made. 

Q II. CIRCULAR-ARC AIRFOILS IN DRY AIR 

In operating the tunnel with dry air, we can assume isentropic flow up to the 
conditions where compression shocks appear. The density distribution in the field 
of observation, determined with the interferometer, can be used to determine the 
corresponding Mach number distribution. 

If the values of the perturbation parameter <p x It H* - 1 as evaluated at the air- 
foil are plotted as a function of wing chord £, we find that the values obtained in 
this tunnel are above the theoretical curve given by Spreiter and Alksne [Ref. 4] 
for an infinitely expanded flow field. The correspondence of these measured values 
with the measurement results from other authors [Refs. 5, 6] is not fully satisfac- 
tory. 

Simultaneously with the optical density measurements, the static pressure was 
measured at the airfoil for 5 ■ 0.25, 0.5, and 0.75. If pressure and density meas- 
urements are compared with each other at corresponding points, then we find syste- 
matic deviations between the measurements, if we assume that the local supersonic 
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region is enclosed by a compression shock sufficiently far down stream. The primary 
reason for these deviations is the boundary layer that exists at the observation win- 
dows. The static pressure in the boundary layer at the tunnel windows remains con- 
stant in a direction normal to the flow and along the direction of the light, while 
the temperature rises to approximately 80 percent of the temperature of the air at 
rest. Based on the equation of state, the density in the boundary layer must there- 
fore decrease. This means that there is a small additional deflection of fringes in 
the interferometer image in the direction of increasing fringe deflection. The inter- 
ference fringe image, when evaluated, leads to higher Mach numbers. 

Thickness and profile of the boundary layer at the window panes are unknown. 
Correction values can be estimated by using the assumption that the boundary layer 
over the tunnel windows is the same in thickness and profile as that which exists 
upstream of the airfoil without boundary -layer suction. 

Applying this correction produces good agreement between pressure and density 
measurements at the pressure-measurement points on the wing. Following correction, /2 00 
the differences are about 2 percent. In the <{> x , 5 diagram, applying this correc- 
tion results in displacing the curves by 45 % 0.03 to the right. Applying the cor- 
rection improves the agreement of the measured values with the results from other 
publications [Refs. 4-6] (Figure 5). 

In order to apply the measured values obtained for the choked tunnel to condi- 
tions of airfoil flow in free space, it is necessary to apply a tunnel correction, 
whose computation is difficult in the transonic region and has so far been possi- 
ble only for special cases [Ref. 7]. Consideration of these special cases indi- 
cates that values measured for a choked tunnel under conditions when the tunnel is 
wide enough (H/L _> 2, t <. 0.1) are in better agreement with values for the flow in 
free space than those obtained with the aid of perforated or slotted transonic tun- 
nels [Ref. 8]. 


5 


The problem of tunnel corrections will not be considered further. The measured 
values obtained here for dry air will only be used for comparison with those for 
humid air and are compared to known measurement results only for control purposes. 

Tests using a circular-arc airfoil glued in mid-tunnel were not carried beyond /201 
the preliminary stages, because the small size of the models necessitated by the 
small tunnel height allowed only the development of a laminar boundary layer over 
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the model. The shock configuration for an airfoil with laminar boundary layer dis- 
plays quite a different picture from an airfoil with turbulent boundary layer (Fig- 
ures 7, 8) [Ref. 9]. Thus, an additional and undesirable parameter would have to 
be considered in the studies. It is for this reason that the studies were only made /204 
with a half-model that was located at the tunnel wall and that always had a devel- 
oped turbulent boundary layer. 

One other phenomenon should be pointed out, which was observed for airfoils 
located ’n the middle of the tunnel. When the incident flow Mach number approaches 
the choke Mach number, nonstationary processes readily appear in the test section. 

They show up as an alternating movement of the shock waves enclosing the local 
supersonic region (Figure 9). Presumably, the tunnel alternately chokes, a phenome- 
non which is initially expected in the vicinity of the test section choking and 
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which is easily triggered through minute asymmetries in the flow field (for exam- 
ple, asymmetrical separation at the airfoil or small angles of attack of the air- 
foil). 


III. CIRCULAR-ARC AIRFOIL IN HUMID AIR 

In contrast to the tests with dry air, the adjustable nozzle located downstream 
from the test section was fully opened and left that way for all tests using humid 
air. The humidity of the air was varied in the range from 0 <_ x < 0.016 kg water/ 
kg air using mixing inside the balloons located upstream from the intake opening. 
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The Lnterferograms for these tests clearly show the changes in density distribution 
in the flow Held for increasing values of humidity (Figures lOa-lOh) . The shock 
wave bounding the local supersonic region moves slowly upstream with increasing 
humidity and becomes weaker* so that* after a certain humidity is reached, the shock 
can no longer reach the opposite wall of the tunnel. Up to this point, the changes 
are quite similar to those observed for dry air for decreasing incident flow Mach 
number. For the tests made with humid air* the incident flow Mach number decreases 
only slightly from M» £ 0.83 for x ■ 0 kg water/kg air down to M w £ 0.81 with x • 
0.01. When humidity increases to reach a value of x £ 0.007, a second shock is 
formed, which does not reach the airfoil at the point upstream from the shock bound- 
ing the supersonic region. The new shock is at times somewhat longer than the one 
that follows it and could be the result of strong condensation phenomena within the 
local supersonic region. Ackeret, Feldmann, and Rott [Ref. 9] observed similar 
double-shock configurations for a curved wall and suspected water vapor condensation 
as the cause. 

When the humidity of the air finally exceeds the value of x ^ 0.012, the shock 
system becomes nonstationary (Figure 11). Shock II, which bounds the local super- 
sonic region, begins to move upstream starting with its whip-like end away from the 
airfoil, while Shock I, which lies further upstream, becomes nonstationary and grows 
weaker as it moves upstream and finally disappears completely. After Shock II has 
moved some distance upstream, it also becomes weaker with the flow accelerating 
behind it, and a new shock, enclosing the local supersonic region, is formed; this 
is Shock 11^. The original Shock II now takes on the role of the roomer Shock I, 
and thus becomes Shock Ii. This process continues as we described above and is per- 
iodic. 
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If one follows the line of constant density p/pg ■ 0.6339, which, for upstream 
isentropic flow, has a Mach number of 1 and thus represents the sonic line for these 
assumptions, we observe the following (Figures 12, j 3) t 

For the range 0 £ x < 0.009, the sonic line thus defined reaches the opposite 
wall at a point slightly downstream. Only downstream from the sonic line can we 
notice changes in the density distribution due to increasing humidity. 

For x % 0.009, the picture changes in that the line p/pg - 0.6339, which starts /206 
from the airfoil at £ ^0.25, no longer reaches the opposite wall. For air humidi- 
ties above x ^ 0.009, it is bent back toward the airfoil in the direction of the 
airfoil wall, after which it again bends and reaches the wall opposite the airfoil. 

This occurs if there is no collision with a shock. The line p/Pq ■ 0.6339 thus 
^sumes an S-shaped form, and we can state with confidence that this line no longer 
r it responds to I! * 1 following the first deflection. The bending of the line 
p/pg ■ 0.6339 is a result of heat addition through condensation downstream from the 
first part of this line. The f]jw in the condensation area is no longer isentropic, 
thus making it possible to determine the sonic line without additional data on heat 

« j, 

addition. |f 

Should the density distribution (p/pg) along the airfoil be plotted as a func- 
tion of £, the following picture evolves: For increasing humidities x, the shock 

bounding the local^j^ersonic region moves upstream (Figure 6). Ahead of the shock, 
all density curves lie within the statistical uncertainties given by the bounds of 
measurement accuracy. Within this range, the effect of condensation is not notice- 
able. Even the double-shock configuration, as shown in Figure lOf, makes itself 
noticed by the strong flattening of the density decrease immediately upstream from 
the shock (Figure 6, curve x - 0.0093). The addition of heat in the field above 
the airfoil has only a minor effect on the density distribution across the airfoil 
upstream from the shock. 
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The family of curves evaluated for humid air strongly resembles that obtained 
for dry air (cf. Figures 5, 6). For the dry air tests, the incident flow Mach num- 
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ber was varied by using the second constriction in the range from M cr itical ^ 
to ^choke % while for the test with humid air, as we have mentioned, the inci- 

dent flow Mach number was changed only due to increasing humidity in the range 

M £ 0.83 to M % 0.81. Thus, the moving ahead of the shock as a function of 
00 00 

increasing humidity cannot be explained by the reduction of the incident flow Mach 
number. The decrease in Mach number upstream from the shock as a consequence of 
heat addition through condensation is probably its chief cause. The effect of 
increasing the amount of heat added that accompanies increasing air humidity is 
clearly shown in the shrinking of the local supersonic region. The fact that the 
incident flow Mach number M^, is barely changed with increasing amounts of heat added 
indicates that the predominant portion of the heat addition occurs downstream from 
the sonic line. Heat addition in the wind tunnel upstream from the sonic line 
results in a decrease of the flow mass. In this case, the incident flow Mach num- 
ber would thus have to decrease. The slight decrease observed in the incident flow 
Mach number for increasing air humidity could be caused by weak condensation phe- 
nomena upstream from the sonic line. The density measurements carried out make it 
impossible to determine the Mach number under anisentropic flow conditions; there- 
fore, we can only guess what happened here. 

Above x % 0.012, the flow becomes nonstationary, especially in the flow field 
above the airfoil (Figure 11). The pressure oscillations caused by the moving 
shock waves are, surprisingly, only barely noticeable at the airfoil itself. The 
airfoil is reached only by the shock bounding the local supersonic region. This 
shock exhibits such small motions at its point of attachment that we can say that 
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Che motions are within the range o£ statistical uncertainties of the measurement 
itself. The density distribution over the airfoil upstream from the shock is also /208 
almost invariant up to a short distance upstream from the shock itself. The non- 
stationary phenomena thus take place only in the flow field and have very little 
effect on the flow conditions at the airfoil itself. This result could not have 
been suspected from the Schlieren photographs, which clearly show a shock comigu- 
ration at variance with photographs made when dry air was used an the medium in the 
wind tunnel (Figures 1, 2). 


IV. CONCLUDING RFMARKS 

I would like to thank Professor Zierep and Dr. Euteneuer for their lively inter- 
est and many suggestions during the tests. 

I also thank the German Research Society for making the interferometer available 
to me along with the optical system and spark photography equipment used with it. 


BIBLIOGRAPHY 

1. Steffen, H.-H. , "Interferometric Density Measurements in Airflows Through Laval 
Nozzles with Energy Addition Using Water Vapor Condensation," Thesis, Technical 
University at Karlsruhe, 1,961. 

2. Schmidt, B. , "Observations of Perturbations Caused by Water Vapor Condensation 
in a Supersonic Windtunnel Nozzle," Dissertation, Technical University at Karls- 
ruhe, 1962, 1962 Yearbook of WGLR, pp. 160-167. 

3. Head, R. M. , "Investigations of Spontaneous Condensation Phenomena," Disserta- 
tion, Calif. Inst, of Techn. , Pasadena, Calif., 1949. 

4. Spreiter, J. R., and Alksne, A. Y. , NACA Rep. 1359, 1958 (supersedes NACA T. N. 
3970). 

5. Knechtel, E. D. , NA FA TN D-15 ( 1959). 

6. Michel, R. , Marchand, F., and LeGallo, J., ONERA Publ. No. 65, 1953, and ONERA 
Publ . No. 72, 1954. 

7. Marschner, B. W. , Journ . Aero Sci . , 23, No. 4, 1956. 

8. Spreiter, J. R. , Smith, D. W. , and Hyett, B. J, , NASA T. R. R-73, 1960. 

9. Ackeret, J., Feldmann, F., and Rott, N., "Research on Compression Shocks and 
Boundary Layers in Rapidly Moving Gases," Communications of the Institute for 
Aerodynamics of the ETH, Zurich, Issue No. 10, 1946. 




ORIGIN AUrAGBB 
OF POOR QUAUXlt 


10 










1 



Figure 3. Schematic Diagram of Wind Tunnel. Dimensions are in meters. 


Key : 

A. Intake tube of 0.175 m diameter from the balloons (18 m^) 

D. Static-pressure measuring-stations 
G. Boundary-layer suction 
Gl. Straightener 
K. Tank (34 m^) 

KH. Tank valve of 0.2 m diameter 
M. Test section 

P Q . Pitot tube for total pressure measurement 
Sb. Screen 

Sch. Gate valve of 0.3 m diameter 

T q . Temperature-sensing element for measuring temperature at rest conditions 
VD. Adjustable nozzle 
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Figure 4. Test Setup as Seen From the Camera Side 


Key : 

1. Windtunnel intake tube 

2. The NVi 200 ball valve 

3. Interferometer and enclosed lest section 

4. Concave mirror 

5. Interferometer control console 

6. Drum Camera 

7. Completed adjustable nozzle 
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Figure 11. [Flow with] T = 10%, H = 148 mm, L = 75 mm, x = 0.017 kg water /kg a • r , 

0.81, Boundary Layer Sucked off Upstream of Airfoil. [Film speed] 
10,000 frames per second 




suction station 

Figure 12. Density distribution for x = 0, t = 10%, H = 148 mm, L = 75 mm, and 
M £ 0.83. 
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Figure 13. Density distribution for X ^ 0.010, t = 10%, H ■ 148 mm, L 
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